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Abstract 

The phase convention independent measure of CP violation for three gener- 
ations of leptons is evaluated at different energy scales. Unlike in the quark 
sector, this quantity does not vary much between the weak and the grand uni- 
fication scales. The behavior of the measure of CP violation in the Standard 
Model is found to be different from that in the extensions of the Standard 
Model. 
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with the discovery of neutrino mixing, which indicates nonzero neutrino masses , there 
recently has been a lot of interest in neutrino physics, including measuring possible CP vio- 
lation in the leptonic sector of the three-generation Standard Model 0. CP nonconservation 
in the leptonic sector has important implication in leptogenesis in the early universe, as it 
can feed into a baryon asymmetry through the sphaleron mechanism [Q. One of the neces- 
sary conditions to have baryogenesis or leptogenesis is CP violation |Q. It is thus natural 
to discuss CP nonconservation in the leptonic sector at higher energy scales. 

In analogy with the quark sector, the Maki-Nakagawa-Sakata (MNS) matrix describes the 
mixing among the flavor eigenstates z/^ and mass eigenstates z/j of neutrinos with a = e, fi,T 
and i = 1, 2, 3 [| 



MNS 



( 

-C23S12 - Ci2S23Si3e*'^ C12C23 - Si2S23Si3e*'' C13S23 
y S12S23 - Ci2C23Sl3e*'' -C12S23 - C23Si2Si3e*'' C23C13 j 



(1) 



where Cjj = cos 6ij, Sij = sin 6ij (0 < 6ij < 7t/2), and 5 is a phase. Instead of referring to the 
above phase-dependent matrix, it has been proven that with three generations of fermions 
CP is not conserved if and only if the single parametrization- independent quantity 

detC = det[NN\EE^ (2) 

is nonvanishing 0, where and E denote the 3x3 Yukawa coupling matrices for the 
neutrinos and charged leptons, respectively^ In terms of Vmns matrix elements, 

detC = 2{ml - ml){ml - ml)[ml - ml)[m\ - mlj{ml^ - mlj{ml^ - ml^)J, (3) 

where m^^^^^^^ and rrir^^^e are eigenvalues of N and E, respectively, and J = 
■Si2'523Si3Ci2C23cf3 sin 5 is an invariant proposed by Jarlskog in the quark sector [||. The 



^As pointed out in Ref. the definition in Eq. (^) implies that the neutrinos are Dirac-type 
particles. If one has Majorana neutrinos, the definition should be detC = det[N^ N, EE^ which, 
however, does not affect later analyses. 



quantity det C vanishes if (i) any two leptons with the same charge have the same mass; (ii) 
any of the angles Oij has the value or 7r/2, or (iii) 5 = or ttJ^ 

To study the scale dependence of det C, it is natural to solve its renormalization group 
equation (RGE) given a set of initial values for the relevant quantities such as the Yukawa 
and gauge couplings at the weak scale. Unlike the situation in the quark sector where the 
analogous quantity drops by 4 to 5 orders of magnitude and may provide insufficient CP 
violation for baryogenesis we will show that this quantity for leptons does not change 
very much between the weak and grand unification theory (GUT) scales, which we take to 
be ~ 10^^ GeV. Since CP violation is one of the ingredients for baryogenesis or leptogenesis 
1^, this may have an important consequence for the evolution of the early universe. 

We will consider two methods for neutrinos to acquire mass: (1) the neutrinos have the 
Dirac masses generated through the the same mechanism as the other fermions; and (2) 
the neutrinos obtain Majorana masses through the seesaw mechanism [0. We will study 
the situation in the Standard Model (SM), the two Higgs doublet model (THDM) and the 
minimal supersymmetric model (MSSM). 

As shown in Ref. |]^ for quarks, det C obeys a rather simple equation under the renormal- 
ization group. In the leptonic sector, this equation can be derived from the renormalization 
group equations for and E. They can be written in the following form: 

^ = [-Gn + Tn + cuNN^ + cuEE^] N, (4a) 

dE r 1 

= [-Ge + Te + C2iEE^ + C22NN^\ E, (4b) 

where t = (l/lGvr^) ln(yu/M^), and Gn, Ge, Tn, Te and Cij are hsted in Table I with U and 
D being the 3x3 up-type and down-type Yukawa coupling matrices, respectively^. Using 
Eqs. (0) and (^), one obtains 



^These conditions are not independent of one another; see, for example, Ref. Q. 

^The trace part of radiative corrections comes from the fermion loops on the Higgs fields where 
quark contributions enter. 
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Models 
Gn 
Ge 
Tn 
Te 

Cll 
Cl2 
C21 
C22 



SM 

20^1 

9n2 , 9 2 
i^l + 4^2 



Tr 
Tr 



St/t/t + ZDD^ + iViVt + EE^ 



Tr 
Tr 



THDM 

, 9 2 
20^1 ^ 4^2 

9„2 _L 9„2 
iSl + 192 

3C/C/t + iViVt 

3 
2 

1 
2 

3 
2 

1 

2 



Tr 
Tr 



MSSM 
3C/C/t + iVM 

3L>L>T + £;£;t 

3 
1 
3 
1 



TABLE I. 



(5) 



where A = -{Gn + Ge) + (T^ + Te) + (2cn + Cu)NW + (2c2i + C22)E^E. Since 

dlndetC 



2TrA 



we obtain a simple evolution equation 

(detC)^ 



(detC) 



exp 



Mz 



(2TrA)dt 



(6) 



Assuming a hierarchy structure in the Yukawa couplings, one can neglect contributions from 
D, N and E as unimportant because at the leading order of the perturbative correction the 
top quark loop in the Higgs field renormalization dominates. Some comments are in order: 

• Assuming that all contributions to the Yukawa coupling evolution are due to a heavy 
top quark, we find that all lepton Yukawa couplings scale in the same way 
Therefore, the ratio {AM"^)^/ {AM'^)mz is independent of the lepton masses, where 
AM2 = {ml - ml){ml - ml){ml - ml){ml^ - mlj{ml^ - mlj{ml^ - ml^). 

• Most of the running of det C is due to the running of masses (instead of that of the 
Jarlskog invariant J), which is mainly a result of the dominant top quark Yukawa 



coupling. By comparing the running of detC and AM^, we find that, as expected, J 
does not change very much over the energy scale range considered. 

When running det C in the three models considered, the vacuum expectation energy of 



the Higgs field is taken to be f = 174 GeV in the SM and v = ^v^ + vj with tan/5 = 
Vu/vd = 10 3 in both the THDM and MSSM. With yt being the Yukawa coupling of the top 



quark, one has [11,12 



-3 {%9l + U) + I8y^ (SM), 
TtA ^ I _3 (l^f + Isl) + 9y^ (THDM), (7) 
-3{fgl + Qgl)+9y^ (MSSM). 

As one can see from Eqs. (|^), the evolution of detC for leptons differs from that for quarks 
because it does not involve the QCD coupling, which significantly drags down the running 
quark masses. 

The evolution of the ratio (det C)^/ (det C) function of the energy scale yU for Dirac 

neutrinos is shown in FIG. 1. It is seen that detC has a very distinct behavior in the SM 
that is different from that in the THDM and MSSM. It is the result of a stronger dependence 
on yt for TrA in the SM in Eq. (|^) and the evolution of yt itself. A convex curve for the SM 
follows from yt becoming smaller at larger scales, thus decreasing TrA. The minor difference 
between the THDM and MSSM is due to the differences between the coefficients of gi and 
(?2- Therefore, the dependence on yt is the dominant factor in determining the evolution 
of (detC)^/(detC)Mz- Within the SM, detC becomes larger as the energy scale goes up; 
whereas in both the THDM and MSSM, det C evolves to a smaller value. At the GUT 
scale, det C in the SM is ~ 6 times bigger than in the THDM and ~ 12 times than in the 
MSSM. The curves for the scale dependence of AM^ are almost identical to those for det C. 
Numerically, (detC, AM^) at the GUT scale is (2.42,2.42), (0.40,0.39), and (0.21,0.19) in 



^The value of tan/3 cannot be too large if the top quark effect is considered to be dominant. If 
tan (3 is very large, bottom quark contributions to the evolution should be taken into account too. 
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FIG. 1. The evolution of (det C)^/(det C)mz with the scale fi. We assume Dirac masses for 
neutrinos. The solid curve is for the Standard Model, the dotted curve is for the two Higgs model, 
and the dashed curve is for the minimal super symmetric model. 

the SM, THDM and MSSM, respectively. 

If a smaller value of tan/? is used, the ratio (det C)^/(det C)m2 iii the THDM and 
MSSM will behave more like that in the SM. This is because in this case yt is larger at the 
electroweak scale compared to the case with a larger tan/?. One should also note that in 
some circumstances the evolution of yt blows up because within the dominant top Yukawa 
coupling assumption, 

V,(r) = ^ "TrL,. (8) 

1 -cyt(O) Jo E{t)dt 

where E{t) = exp [— Jq ai gfdt] and c > 0. The numerical values of a, with i = 1, 2, 3 and 
c depend upon the model. Eq. tells us that either one has to include other fermion 



contributions or that higher order perturbative corrections ||T3|] should be considered when 
the denominator gets close to zero. 

The other interesting case occurs when the neutrino masses are acquired through the 
seesaw mechanism. In this case, we have the effective interaction Nip^^ tp^^ H H. This will 
introduce two extra loop contributions on the Higgs legs compared to the previous case. 
Therefore, in the seesaw model both Gn and T/v are twice as big as for Dirac neutrinos. 
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FIG. 2. The evolution of (det C)^/(det C)mz with the scale /i. The labeling convention is the 
same as FIG. 1, but neutrino masses are generated through the See-Saw mechanism. 

We see in Fig. 2 that the quantity det C at the GUT scale is roughly twice as big as in 
Fig. 1. (detC,AM2) at the GUT scale is (5.05,5.05), (0.84,0.79), and (0.49,0.42) in the 
SM, THDM and MSSM, respectively. 

Regardless of their different behavior, we see that for all models considered the evolution 
of det C is, unlike the quark sector 0, of order 10~^ to 10 in going from the electroweak scale 
to the GUT scale. The actual factors that enter a CP-violating quantity are not necessary 
only det C, but may involve other mass scales relevant to the process. The CP-violating 
mechanism responsible for leptogenesis may have a different origin. As long as the physics 
the gives rise to lepton asymmetry effects at a higher energy scale is in the same multiplet of 
the Standard Model particles within the unified theory, as an effective theory, the quantity 
det C will still be a good indicator of how much CP nonconservation would be at higher 
energy scales, though the exact relation depends upon models. 



As delineated in Ref. fl^ with a toy model containing two very heavy SU(2)-doublet 
scalars, if the neutrinos acquire small Dirac masses, the CP- violating processes due to 
the heavy scalars at the large energy scale can produce a neutrino asymmetry. Both the 
heavy scalars and the light Dirac neutrinos are out of equilibrium at the weak scale where 



sphalerons come into play to produce baryon and lepton number violation. Similarly, if 
the neutrino masses are generated through the seesaw mechanism with very heavy (~ 10^^ 
GeV) singlet right-handed neutrinos, then both the out-of-equilibrium and CP violation 
conditions can be satisfied around the GUT scale. One then can have lepton asymmetry 
leading to baryogenesis [Q. 

We conclude that hierarchical neutrino masses suggest that the possible leptonic CP 
violation strength at the GUT scale differs from that at the electroweak scale by no more 
than one order of magnitude. Since this quantity is reparametrization invariant, models of 
grand unified theories should give consistent predictions for it. The fact that the measure of 
the leptonic CP nonconservation does not have a significant decrease at higher scales may 
help leptogenesis in the early Universe. 
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of Energy under Grant No. DE-FG02-91ER40682. 



8 



REFERENCES 



[1] Super-Kamiokande Collaboration, Y. Fukuda et. al, Phys. Rev. Lett. 81, 1562 (1998); 
Super-Kamiokande Collaboration, K. Scholberg, e-Print: piep-ex/9905016| (1999); F. 
Ronga, Nucl. Phys. Proc. Suppl. 87, 135 (2000). 

[2] V. Barger, K. Whisnant, and R.J.N. Phillips, Phys. Rev. Lett. 45, 2084 (1980); H. 
Minakata and H. Nunokawa, Phys. Lett. B 413, 369 (1997); M. Tanimoto, Phys. Rev. 
D 55, 322 (1997); K. Dick, M. Freund, M. Lindner, and A. Romanino, Nucl. Phys. B 
562, 29 (1999); M. Koike and Joe Sato, Phys. Rev. D 61, 073012 (2000); Erratum-ihid. 
Phys. Rev. D 62, 079903 (2000); V. Barger, S. Geer, R. Raja, and K. Whisnant, Phys. 
Rev. D 62, 073002 (2000). 

[3] M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986); M.A. Luty, Phys. Rev. D 
45, 455 (1992); E. Ma and U. Sarkar, Phys. Rev. Lett. 80, 5716 (1998); R. Barbieri, P. 
Creminelli, A. Strumia, and N. Tetradis, Nucl. Phys. B 575, 61 (2000). 

[4] A. Sekharov, Pisma Zh. Eksp. Theor. Fiz. 5, 32 (1967); JETP Lett. 5, 24 (1967). 

[5] This mixing matrix can be parametrized in various ways; see, for example, Z. Maki, 
M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962); J.N. Bahcall, M. 
Kamionkowski, and A. Sirlin, Phys. Rev. D 51, 6146 (1995); J.N. Bahcall and P. 
Krastev, Phys. Lett. B 436, 243 (1998); R. Escribano, J.M. Frere, A. Gevaert, and 
D. Monderen, Phys. Lett. B 444, 387 (1998). 

[6] C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985); Phys. Rev. Lett. 57, 1875 (1986); Z. 
Phys. C 29, 491 (1985); Phys. Rev. D 35, 1685 (1987). 

[7] Z.Z. Xing, e-Print: ^ep-ph/0009294| (2000). 

[8] H. Fritzsch and Z.Z. Xing, Nucl. Phys. B 556, 49 (1999). 

[9] G.G. Athanasiu, S. Dimopoulos, and F.J. Gilman, Phys. Rev. Lett. 57, 1982 (1986). 

9 



[10] T. Yanagida, in Proc. of the Workshop on the Unified Theory and Baryon Number in the 
Universe (KEK report 79-18), edited by O.Sawada and A. Sugamoto, p. 95 (1979); M. 
Gell-Mann, P. Ramond, and R. Slansky, in Supergravity, edited by P. van Nieuwenhuizen 
and D.Z. Preedman, p. 315, North Holland, Amsterdam (1979). 

[11] B. Grzadkowski and M. Lindner, Phys. Lett. B 193, 71 (1987); B. Grzadkowski, M. 
Lindner, and S. Theisen, Phys. Lett. B 198, 64 (1987). 

[12] T.P. Cheng, E. Eichten, and L.F. Li, Phys. Rev. D 9, 2259 (1974); V. Barger, M.S. 
Berger, and P. Ohmann, Phys. Rev. D 47, 1093 (1993); H. Fusaoka and Y. Koide, 
Phys. Rev. D 57, 3986 (1998); P. Kielanowski, S.R. Juarez W., J.G. Mora H., Phys. 
Lett. B 479, 181 (2000). 

[13] D.R.T. Jones and L. Mezincescu, Phys. Lett. B 136, 242 (1984); P. West, Phys. Lett. B 
137, 371 (1984); A. Parkes and P. West, Phys. Lett. B 138, 99 (1984); J.E. Bjorkman 
and D.R.T. Jones, Nucl. Phys. B 259, 533 (1985); V. Barger, M.S. Berger, and P. 
Ohmann, Phys. Rev. D 47, 1093 (1993); H. Fusaoka and Y. Koide, Phys. Rev. D 57, 
3986 (1998). 

[14] K. Dick, M. Lindner, M. Ratz, and D. Wright, Phys. Rev. Lett. 84, 4039 (2000). 



10 



